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Abstract 

We study the lepton flavor models with the flavor symmetry (Zjy x Zn x Zn) x Z3. Our 
models predict non- vanishing discrete values of #13 as well as #12 and #23 depending on N. 
For certain values, our models realize the tri-bimaximal mixing angles with #13 = 0. For 
other values, our models provide with discrete deviation from the tri-bimaximal mixing 
angles. 



1 Introduction 



It is the important issue to study the origin of the flavor structure in particle physics, that 
is, the hierarchy of quark/lepton masses and their mixing angles. The form of lepton mixing 
angles is quite different from one of quark mixing angles (see e.g. for the lepton mixing [IHS]). 
The lepton mixing angles are large, while quark mixing angles are smaller. The tri-bimaximal 
mixing Ansatz in the lepton sector is quite interesting [I] and the exact tri-bimaximal mixing 
leads to 9 13 = 0. 

Non-Abelian discrete flavor symmetry is a powerful tool to derive masses and mixing 
angles of quarks and leptons (see for review Refs. [5JE] and references therein). Especially, 
large mixing angles of leptons can naturally arise by assuming non-Abelian discrete flavor 
symmetries. Indeed, there are many models based on flavor symmetries to explain fermion 
masses and mixing. 

Recently, the results of T2K and Double Chooz reported a non- vanishing value of 9 is [7f8] . 
Thus, it is important to derive the lepton mixing with #13 7^ 0. However, the tri-bimaximal 
mixing matrix is still good at a certain level. Thus, it is interesting to study a small deviation, 
e.g. from the tri-bimaximal mixing angles with 9 r3 7^ 0. After the result of T2K is released, 
several attempts have been studied to to explain non- zero with various methods [9H3T]. 
For instance, introducing somewhat larger corrections to the model of tri-bimaximal mixing, 
one can obtain non- zero value for U e3 . However, such a method contains a free parameter 
so that it does not predict a specific value. In this paper, we build a model that predicts 
a special mixing pattern with deviations from the tri-bimaximal mixing angles by discrete 
values, which are consistent with experiments. 

A simple way to derive the tri-bimaximal mixing is to assume residual symmetries for 
lepton mass matrices. Precisely, one can realize the tri-bimaximal mixing when there is the 
Z 3 symmetry for the charged lepton mass matrix and the Z2 symmetry for the neutrino mass 
matrix. These Z 3 and Z 2 symmetries may be originated from larger non-Abelian flavor sym- 
metries like A4, S4, etc. Thus, it is interesting to use non-Abelian discrete flavor symmetries 
such as A4 and S4 in order to derive the (exact) tri-bimaximal mixing [5j[6] . To predict non- 
vanishing #13, we propose a modified symmetry for the charged leptons. In our approach, 
we use the (Z N x Z N x Z N ) x Z 3 symmetry. With a choice of proper vacuum alignment, 
the charged lepton sector retains Z 3 and Zn symmetries while the neutrino sector has Z2 
symmetry, which is the subgroup of Zn- Due to the additional Zn symmetry of charged 
lepton, we obtain non- vanishing #13. 

This paper is organized as follows. In section 2, we study about group-theoretical aspects 
of (Z N x Z N x Z N ) x Z 3 by giving conjugacy classes, character table, and tensor products. 
In section 3 we introduce the model based on {Zn x Zn x Zn) x Z 3 lepton flavor symmetry 
with one singlet and two triplet flavon fields. In section 4 we ensure the vacuum alignment 
assumed in the model by two ways, that is, analyzing scalar potential and using boundary 
conditions in extra dimensions. In section 5, higher dimensional terms are taken into account. 
Section 6 is devoted to the summary. 
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Z N x Z N x Z N ) x Z x 



We use the group (Zn x Zjv x Zjv) x Z 3 for the lepton flavor symmetry. Here, we study 
group-theoretical aspects of our group. We denote the first, second and third Z N generators 
by a, a' and a" , respectively, and the Z 3 generator is denoted by b. They satisfy the following 
algebraic relations: 



N _ „lN _ „ilN _ u3 _ „ „„'_„'„ „„"_„"„ J Ji _ Jl J 

(1) 

> 



a N = a' N = a" N = b 3 = e, aa' = a'a, aa" = a" a, a' a" = a" a 
ba n b~ l = a //n , ba^b" 1 = a n , ba" n b~ l = a' n 



where e denotes the identity. Group elements can be written by b k a e a >m a" n with £,m,n = 
0, N — 1 and k — 0, 1, 2. When iV/3 = integer, this group is isomorphic to S(3iV 3 ) [B]. For 
-/V/3 7^ integer, this group corresponds to A(3iV 2 ) x Z N . Here, we use both types of groups. 

First, let us study the conjugacy classes of the elements without including b, i.e. a e a' m a" n . 
These elements elements a e a' m a" n satisfy the following relations: 

ba t a fm a" n b- 1 = a m a ,n a" 1 , b 2 a l a' m a" n b- 2 = a n a' e a" m . (2) 

Note that ba^a'^a'^b" 1 = b 2 a e a' e a" e b~ 2 = a e a' e a" e . Thus, we find the following conjugacy 
classes of a i a' m a" n : 

d : {e}, 

Cf : {a W}, for £ = 1, ■ ■ ■ , N - 1, (3) 

Cf m ' n) : {aV"V"\ a m a' n a" e , a n a' e a" m }, for £, m, n = 0, • • • , N - 1. 

The last classes c^' m ' n ^ exclude the case with I = m = n. The number of the classes c^' m ' n ^ 
is equal to (iV 3 — N)/3, while the number of the classes C± is (N — 1). 

Next, we study the conjugacy classes of the elements including b, that is ba l a' m a" n and 
b 2 a e a' m a" n . It is found that 

a p a' q a" r {ba l a' m a" n )a- p a'- q a"- r = ba e+r - p a' m+p - q a" n+q - r , 
ba p a' q a" r {ba e a' m a" n )a- p a'- q a"- r b- 1 = ba m+p - q a' n+q - r a" e+r - p , (4) 
b 2 a p a' q a" r {ba l a' m a" n )a- p a'- q a"~ r b~ 2 = ba n+q ~ r a' i+r - p a" m+p - q . 

Note that the sum of the powers of a, a', and a", i + m + n, does not change. Then the 
conjugacy classes of ba e a' m a" n are obtained as 

C$ : {ba e a' m a" p - e - m \£, m = 0, • • ■ , N - 1}, for p = 0, • ■ ■ , N - 1. (5) 

Similarly, we obtain the conjugacy classes of b 2 a e a' m a" n as 

C'$ : {&Va'"V /p '-^ m |£, m = 0, • • • , N - 1}, for p = 0, • ■ ■ , N - 1. (6) 

Then, the total number of conjugacy classes is 

1 + (N - 1) + (iV 3 - N)/3 + iV + iV = -iV(iV 2 + 8). (7) 
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Now, let us study representations of our group. Suppose that there are m n n-dimensional 
irreducible representations, where group elements are represented by (n x n) matrices. The 
identity e is always represented by the (n x n) identity matrix and its character x( e ) is 
obtained as x( e ) = n on the n-dimensional representation. The orthogonality condition of 
characters requires J2 n n 2 m n to be equal to the order of the group i.e., 

mi + 4m 2 + 9m 3 H = 3iV 3 . (8) 

In addition, the total number of irreducible representations should be equal to the total 
number of the conjugacy classes, i.e., 

m 1 + m 2 + m 3 + --- = -N(N 2 + 8). (9) 

The solution of Eqs. © and <Q is obtained as (m u m 2 ,m 3 ) = (3N,0,N(N 2 - l)/3). That 
implies that there are 3N singlets and N(N 2 — l)/3 triplets. 

It is straightforward to derive the 3N singlet representations, if with k = 0, • • • , N—l and 
£ = 0, 1, 2. Their characters for a, a', a" and b are obtained as Xi e k ( a ) — Xi e k ( a ') — Xi e ( a ") = P k 
and Xi e (b) = where p = e 2m ^ N , and to = e 2m ^ 3 . These are shown in Tabled! where h 
denotes the order of the element g in the conjugacy class, i.e., g h = e. Note that these 
characters are nothing but representations for a, a', a" and b on l k . 

On the other hand, we can write a, a', a", and b as 







a = I (t | , a"=\ p n | , b = \ 1 | , (10) 

p e 

on the triplet 3[^[ m ][ n ], where the notation [£][m][n] means 

[^][m][n] = (£, m, n), (m,n,£), or(n,£,m). (11) 

Their characters are shown Table [TJ 
We write components of triplets as 

( X£ \ 

3[f][m][ra] = Um • (12) 

Subscripts of the components describe charges. The tensor product between two triplets 
is given by 





3 [f'][m'][n'] 




°[e+l'][m+m.'][n + n'] v ' ^ [m+n'] [n+l'\ [I + m>] v ' ' s [n+m'\[t + n'\[m+t'\ 




(13) 



1 See e.g. [B] and references therein. 
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ft 


Xl l 


^ 3 [f][m][n] 


Ci 


1 


1 


3 




TV 

gcd(7V,p) 


ppm 


^p(£+m+n) 




TV 

gcd(iV,p,g,r) 


p{p+q+r)m 


ppt+qn+m.r _|_ ppm+ql+rn _|_ ppn+qm+rl 


/nr(p) 


3N 
gcd(7V,p) 


u e ff m 





✓VW 

iV 2 


gcd(7V,p) 








Table 1: Characters. 

If all the subscripts of 3[ e+e >][ m+m '][ n+n >], 3[ m+Tl /]fc-f-^][*+m'] or 3[ n+m /][^ +ri '][ m+ ^] are the same, 
such a triplet can be decomposed into singlets as 

(x a , x b , x c ) 3[k][k][k] = (x a + x b + x e )io + (x a + u 2 x b + UJX C ) 1 1 + (x a + UX b + W 2 X C )!2. (14) 

The tensor product between singlet and triplet is obtained as 

x (y) x v = I u k 'x m y | (15) 
Ku^'xny, , 

d M[m][n] 7 3[<+fe][m+fc][n+fc] 

The tensor product between singlets is obtained as 

0)i| x = (^IjW • (16) 

R v k+y 




3 Lepton mass and mixing 

We use the group (Z N x Z N x Z N ) x Z 3 as the lepton flavor symmetry. We consider the 
case with N/4 = integer. We also assume the additional Z4 flavor symmetry. We assign 
the quantum numbers of the three families of lepton doublets {(-eA^Ar) and lepton singlets 
(e c , /i c , r c ) and the Higgs fields H u ^ under the flavor symmetry as shown in Table [2j These 
denote chiral superfields. In addition, we introduce the flavon fields, Xoi {Xii X^ X*s) an d 
(xi, xli xi)j which also denote chiral superfields. Their quantum numbers under [Z^ x Z^ x 
Z^) x Z 3 and Z 4 are also shown in Table |2J 





e c fi c t c 


(•^ej ^/i) ^r) 




Xo 






(Z 3 N ) x Z 3 
Z 4 


-il 12 10 

A 3A r /4+p ± 3N/4+q x 3JV/4 

1 1 1 


3[3AT/4][AT/4][JV/4] 

3 


1° 



1° 

± N/2 

2 


3[A72][0][0] 

2 


3[jv-i][o][o] 




Table 2: Quantum numbers of leptons, Higgs fields and flavon fields under (Zjy x Z^ x Z^) x Z 3 
and Z4. Parameters p and g of (Z^r x Z^r x Z^r) x Z 3 charges are arbitrary integers. 



The terms relevant to charged lepton masses in the superpotential are obtained as 

w e =y e e%£ e ( X {) N/2 + Uxi) N/2 + £M) N/2 )HMxki) P /^ N/2+3p 

+ y>Ve(x[) N/2 + + £Axi) N/2 )H d (xlxlxi)VA N/2+3q (17) 

+ yVWxi)^ + Uxi) N/2 + ir(xi) N/2 )H d /A N ^, 
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where A is the cut-off scale and y e ,y' 1 ,y T are dimensionless couplings. Similarly, the terms 
relevant to neutrino masses in the superpotential are obtained as 

+ y U 2 (le^X3+WrX2+WrXl)H u H u /A 2 , [ ' 

where y v l2 are dimensionless couplings. 

We assume that all scalar fields develop VEVs which are written by 

(H Uid ) = v u4 , (x v ) = a„A, ((rf, X2, Xs)) = «, a%)A, ({xi xi xl)) = (<4 <4, oJ)A.(19) 
Furthermore, as the vacuum alignment, we assume 

((Xi, Xa, Xs)> = 0, 0)A, ((xi , X2, Xs)) = ^(1, p, p')A, (20) 

where p and p' are phases of Zjv so that we write p = e 2mm / N ; p' = e 2™/7V w ^jj integers m 
and n. 

Taking the vacuum alignment, mass matrices of charged leptons and neutrinos, M e and 
M u , are given by 

f y e af(pp') p 

M e = v d af /2 ( y^ippy 





'y\oi v 




(21) 



</X ^a, 



Except additional phases p and p', the mass matrices are the same as those of the A± model 
|38j . They can be easily diagonalized. The masses of charged lepton are 

m e = y a/ y v d , m M = y^a t v d , m ll = y a/ v d . (22) 

Taking p « 2g, we can realize the mass hierarchy of charged leptons when the Yukawa 
couplings, y e ,y' 1 and y T , are of the same order each other. Suppose that Yukawa couplings 
are of 0(1). Then, the value of 3p should be taken as the same magnitude of N. 
The mixing matrix becomes 

1 
l/y/2 -l/y/2 
l/y/2 l/y/2 

1 u 2 {p + u 2 p')/y/2pp' u 2 {p-u 2 p')/V2pp'\ 
1 u(p + cup')/V2pp' u(p-up')/V2pp' . 
1 (p + p>)/V2ppi {p-p>)/y/2~pp> J 

Assuming p ps u 2 p', the (1, 3) element can be small enough to be fitted to experiments. Then 
we obtain 

v 2 v 2 v 2 

\mi\ = \v\ol v + y v 2 oi v \^ |m 2 | = \y\a v \^ \m 3 \ = \y\a v - y v 2 ot v \^-. (24) 





(23) 



To explain the neutrino mass difference, we need 



2\yla v U^\{cos A)|| = -^AmL, bXI 2 ^ = - Am 2 ol , (25) 

where A is the phase difference between y\a v and y% a' u . Since experiments show | Am 2 tm |/2 > 
Am 2 ol , the second equation gives Am 2 tm > 0, i,e, the normal mass hierarchy. Using them, 
we obtain the neutrino mass spectrum as 



AmL . 2 ■ Am 2 

i S0l' 



1 11 8(l-2r)cos 2 A so1 ' 1 2| ~ 8(1 - 2r) cos 2 A 



Am 2 (26) 
|m3|2= 8(l-2r)7os 2 A +Am -- Am ^ 

where r = Am 2 ol /Am 2 tm . Inserting the best-fit values for Am 2 tm and Am 2 tm and cos 2 A = 1, 
we obtain \m\\ = 16meV, \ni2\ = 18meV, and |m3| = 52meV. The smallest value of sum of 
neutrino masses is 86meV. 

In the usual notation, mixing angles can be expressed by 



sin 2 # 
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2 + cos(27r(m — n )/N) + y/3 sin(2vr(m - n) /N) 



sin 2 ^ 



6 
2 

12 



4 - cos(27r(m - n) /N) - \/3sm(27r(m - n) /N) ' 

2 4sin 2 (7r(m — n)/N) 

sin c'23 = 1 — 1= . 

4 - cos(27r(m - n)/N) - V3sin(27r(m - n)/N) 

The Dirac CP phase is always vanishing. When N = 8, the smallest |i7 e 3| is obtained as 



(27) 



\U*\ = y(2- V2 + V3)/6, (28) 

which is about 0.107. In this case, we finci 

• a 2 ■ a V5-3V2-V3 + VQ 
sin #i2 = — , sin #23 = — . (29) 

V8 + V2 + V6 + 

For N = 12, U e 3 can be vanishing and mixing matrix becomes tri-bimaximal mixing. 

We have calculated possible mixing angles which can be consistent in experiments up to 
N = 40 in Fig. 1. By fixing p and p', we can predict mixing angles. In Fig. 2, we estimate the 
effective mass of double beta decay m ee = [rriiU^ + m 2 t/ 2 2 + m 3 U^ 3 \. Due to Majorana phases 
and uncertainty of neutrino mass differences, the effective mass is predicted in continuous 
line with some fluctuation. For N = 8, the minimum value of the mass of double beta decay 
is predicted about 5.3 meV, and for other N it is about 4.6 meV. 



1 These mixing pattern is obtained by another approach 
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■0 30 40 °- 000 °- 005 °- 010 0015 0020 0025 0030 

N sin 2 *),, 



Figure 1: The left figure shows the relation between N and \U e ^\ for possible values of p and 
p'. The right figure shows the values of sin 2 12 , sin 2 # 13 , and sin 2 #23 for the region < N < 40 
in the case N/4 is integer. Continuous line is the correlation of mixing angles in the limit of 
N -> 00. 



1.000 




0.02 0.05 0.10 0.20 0.50 10 15 20 25 30 35 40 

m„ lin [eV] N 

Figure 2: The left figure shows the relation between minimum neutrino mass and effective 
mass of double beta decay for N = 8. The right figure shows the minimum value of the 
effective mass of double beta decay for each N. 



4 Vacuum alignment 

Here, we discuss the vacuum alignment of flavon VEVs. We study it in two ways, one is to 
use the potential analysis in four- dimensional field theory and the other is to use boundary 
conditions in extra dimensions [39ti42] . 



4.1 Potential analysis 

Here, we study the vacuum alignment of flavon VEVs by using four- dimensional potential 
analysis. The superpotential terms of flavon fields xt i} = 0, 1,2,3) up to four dimensional 
operators are 

«>/ =M 1 (X U ) 2 + M 2 ((vn 2 + ( X u 2 ) 2 + (Xtf) + A!(Xo) 4 + Mx u mx1) 2 + {xlf + {xlf) 

+ A 3 ((xD 2 + u{x v 2 f + u\xl?)({x\) 2 + uj\x2? + "(Xs) 2 ) , . 

+ H(x\) 2 + {xlf + Od0 2 )((vn 2 + (X2) 2 + (x!0 2 ) 1 } 

+ H(X2Xs) 2 + (XsXi) 2 + (X1X2) 2 ) + A 6 ((xr) 4 + (X2) 4 + (X3) 4 ), 
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where Aj = Aj/A (i = 1, 2, • • • , 6) with dimensionless coupling constants Aj. The conditions 
of the potential minimum are written as 



2M lX & + 4A 1 (^) 3 + 2A 2 (x£)((^) 2 + (Xi? + (Xsf) = 0, 
2M 2X r + 2X 2 ( X U ) 2 X U 1 + 2A 3 xr(2(Xi) 2 - (Xs) 2 - (Xs) 2 ) 
+4A 4X r((xr) 2 + (XsT + (xlf + 2A 5X r((x^) 2 + (Xaf) + 4A 6 (xD 3 = 0, 



+4A 4 X2((XD 2 + (x 2 f + (Xaf) + 2\,X 2 ((x\f + (Xs) 2 ) + 4A 6 ( X 2) 3 = 0, 

2M 2X l + 2A 2 (^) 2 X 3 ' + 2A 3 Xs(-(xr) 2 ~ (x^) 2 + 2{xlf) 
+4A 4 Xs((xD 2 + (X2) 2 + (Xs) 2 ) + 2A 5X ^((xr) 2 + (X2) 2 ) + 4A 6 (xD 3 = 0. 



There is a solution we have used as the vacuum alignment 



The leading and next leading order terms of xf (i — 1,2, 3) in the superpotential are 
obtained as 



w f =m'((x[f + (xl) N + (xl) N ) 

+ K(( X {f + coixif + ^(xlf)((x{f + u?toZ) N + wbi) N ) 

+ K((xif + (xif + (xlf)((xif + (xif + (xif) 

+ K((xixif + (xixif + (xixlf) + K((xi) 2N + (xi) 2N + (xlD, 



where m' = m'/A N 3 and A^ = A^/A 2Ar 3 (i = 1,2,3,4) with dimensionless constants vn! and 
A^. When x{ 7^ 0, x 2 0, an d X3 0, we have 



m' + A' 1 (2( X 7V - (%X " (Xif) + 2A' 2 (( X ^ + (xif + (xif) 

+K((xlf + (xif) + 2K(x{f = o, 

m' + K(-(x[) N + 2(xT - (Xif) + 2A' 2 (( X ^ + (xS" + (xlf) 

+K((xif + (x{f) + 2K(xif = o, 

m' + K(-(xif ~ (Xlf + 2(x^) + 2X' 2 (( X {f + (xif + (xif) 

+K((xlf + (xif) + 2\> 4 (xif = o. 



Moreover, when 3X[ — A3 + 2A4 7^ 0, we have the relations 



Then we can realize the vacuum alignment xi — XiP an d X3 = XiP' ■ 



2M 2X U 2 + 2\ 2 (xl)\ v 2 + 2\ 3X2 (-(xl) 2 + 2(x^) 2 - (Xs) 2 ) 



(31) 




(32) 
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There are the cross terms between x u an d X 1 ■ F° r instance, if (xf/A) 7 ^ 2 ~ Xo ~ Xi> the 
most effective terms are 

HKixoT + K((xiT + (x!D 2 + (x^) 2 ))((x £ 1 ) 7V + (xl) N + 

+ KM(xlxt) N/2 + X^X^ 2 + Xa^X^ 2 ) (36) 

+ K(x2Xs(xlxt) N/2 + XixMxt) N/2 + XiX2(xtxt) N/2 ). 

These terms would violate the vacuum alignment that we obtain above. We need to assume 
these cross terms should be suppressed sufficiently in this system. 



4.2 Vacuum alignment on orbifold 

Here, we discuss another way to realize the vacuum alignment for (xi-> X21 X3) an d (Xi> X25 X3); 
using extra dimensional field theory. We study the flavor symmetry breaking by boundary 
conditions on the orbifold [39TH2] M 

We consider eight- dimensional field theory on the T 2 jZ 2 x T 2 /Z 3 orbifold. The fifth and 
sixth dimensions, (x 5 ,x 6 ), are compactified on the T 2 /Z 2 orbifold and the seventh and eight 
dimensions, (x 7 ,x 8 ), are compactified on the T 2 /Z 3 orbifold. It is useful to denote the extra- 
dimensional coordinates (x 5 , x 6 , x 7 , x 8 ) by a complex space z = x 5 + ix 6 and z' = x 7 + ix 8 . 
We assume that (Xi; X25 X3) &nd (xi> X25 X3) ii ve on T 2 /Z 2 and T 2 /Z 3 , respectively. We have 
to fix the boundary conditions of these fields under the Z 2 twist P and the Z 3 twist R, i.e. 

'xKcuz'Y 

X i(uz>) I = P I X2V) I , (37) 





yXl^Z 1 ) 




where u = e 27 ™/ 3 . 
When we take 



-1 



K ■ 

only the component Xi nas the zero-mode. Furthermore, when 

/0 // 0\ 

i? = p m , (39) 
\p n / 

where £ + m + n = mod iV, only the direction (xiiXsLxIj — Xi(l> P* \ P*)i where i,j = 
0, 1, • • • , N—l, has the zero- mode. If these zero modes develop their VEVs in four-dimensional 
effective field theory, the vacuum alignment can be realized. 



The orbifold models are also interesting as the origin of non-Abclian discrete flavor symmetries 
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5 Higher order corrections 

Let us consider higher dimensional terms for charged leptons. We do not consider terms with 
(Xo) 2 or (xi, X25 X3) 2 because it just changes the mass eigenvalues. Next-to-leading terms 
which modify eigenvectors are 

Aw e = e c (4, £», £ T )H d (Ay e lX "o(Xi, & Xs) + ^(Xi, X2, X u 3 ) 2 Mxlx^) N/2+a /^ N/2+3a+2 

+^(l e J tM J T )H d (Ay»Xo(xlx U 2,X^) + >& xtfMxixl) N/2+ V^ N/2+3b+2 (40) 

+r c (4, £ M , ^^(AyrxSOti, X5, Xa) + Ay 2 r (;tf , *5) WiX^WA w/a+a . 

Compared to the leading order, corrections are of 0{afa1). 
Similarly, correction terms for neutrinos are 

Aw„ =(4, l„ 4)(4, € MI 4)(A^(xr, X2, Xs) + A^)(xks, X1X3, X {xi) N/2 H u H u /A N+2 . (41) 

Corrections are of 0((a^) ) to the leading terms. This can be estimated as (m T /y T Vd) 2 - If 
tan /3 is not large, the corrections can be suppressed. 



6 Conclusion 

We have studied the models with the lepton flavor symmetry [Z^ x x Z^) x Z 3 in order 
to extend the tri-bimaximal mixing. The tri-bimaximal mixing can be derived by residual 
Z 3 and Z 2 symmetries for charged lepton and neutrino mass matrices, respectively. We 
extend the symmetry of charged leptons to Z3 and Zn symmetries, then non-vanishing 6*13 
can be obtained. For example, when we choose N = 8, our model predicts sin #13 « 0.107, 
sin6> 12 « 0.337, sin #23 ~ 0.651, and <5cp = 0. For larger values of N, we predict other sets of 
mixing angles, but in any case, our model predicts 5qp = 0. 

We have also predicted the mass spectrum of neutrinos. Three neutrino masses are 
expressed by two complex parameters, so that they are correlated. In our model, only normal 
mass hierarchy is allowed. To compare with future experiments, the relation of effective mass 
of double beta decay and minimum neutrino mass are calculated. 
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